typeset using JPSJ.sty <ver.l.Ob> 



o 
o 
o 

(N 

CD 
00 

0^ 






Magnetization plateaus and phase diagram in polymerized S = 1/2 

XXZ chains 



3nQ, 



\^*\ 



Wei Chemj, Kazuo HiDA^ta and B. C. Sanctuary 

Department of Chemistry, McGill University, Montreal, PQ, Canada H3A 2K6 
^Department of Physics, Saitama University, Urawa, Saitama, Japan 338-8570 

(Received February 1, 2008) 

The magnetization plateaus of p-merized S — 1/2 XXZ chains are studied for general values of 
p. Two plateau-non-plateau critical lines and one plateau-plateau critical line are found for each 
value of p. The universality class of the plateau-non-plateau transition belongs to Brezinskii- 
Kosterlitz-Thouless (BKT) type and that of the plateau-plateau transition, to the Gaussian type. 
The critical points are determined by level spectroscopic analysis of the numerical diagonalization 
results for 4 < p < 8. The multicritical points are calculated using the integral equations based on 
the Bethe ansatz solution of the XXZ model. The behavior of multicritical points are analyzed in 
detail for large p. It is found that the plateau region is enhanced with the increase of periodicity 
p although the non-plateau region persists as far as p is finite. 

KEYWORDS: magnetization plateau, level spectroscopy, exact diagonalization, BKT transition, twisted boundary 
condition 
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§1. 



Introduction 



One-dimensional quantum spin systems are j-attract- 
ing both experimental and theoretical interest .EP Vari- 
ous phenomena have been discovered which result from 
strong spin-spin correlations and strong quantum fluc- 
tuations. Among these phenomena, the magnetization 
plateaus are attracting broad interest as an essentially 
macroscopic quantum phenomenon in which macroscopic 
magnetization is quantized to fractional values of the sat- 
urated magnetization value. Such plateaus are predicted 

and /or pbsprvfjH JTi rnanv Inw dirriensinnal spin SyS- 

tems.BBBii3l'BBBBM00S) The plateau state 
can be also regarded as a spin gap state with non- 
zero magnetization. Oshikawa, Yamanaka and AffleckEP 
proposed the necessary condition for the magnetization 
plateaus as p{S — m^) — q = integer, where p is the 
periodicity of the magnetic ground state in the thermo- 
dynamic limit, S is the magnitude of the spin and ?7i^ is 
the magnetization per site. 

As seen from the above criterion, the magnetization 
plateau appears due to the periodic superstructure of the 
ground state. Therefore it is of interest to investigate 
how the plateau develops with the periodicity p for a 
specific model. In the present work, we investigate this 
problem concentrating on the highest plateau {q — 1) 
of one dimensional p-merized S = 1/2 XXZ model for 
general values of p. 

This paper is organized as follows. In the next sec- 
tion, the model Hamiltonian is defined and the numeri- 
cal results of the phase diagram is presented. For each 
value of p, we find one Gaussian plateau-plateau critical 
line and two plateau-non-plateau Brezinskii-Kosterlitz- 
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Thouless (BKT) critical lines. Three critical lines merge 
at a multicritical point. The p-dependence of the mul- 
ticriticaLncuiiL is investigated in detail by Bethe ansatz 
methodB'B'B in § 3. The final section is devoted to a 
summary and discussion. 

§2. Numerical Results 

2.1 Model Hamiltonian 

The Hamiltonian of p-merized 5=1/2 XXZ chain in 
the magnetic field is given by 

L p-2 
1 = 1 i=0 

+ [i + {p- i)t]npi+p-i,pi+p{A)} (2.1) 

pL 



5MBi/$^5f, 



1=1 



where 

Wm+i(A) 



QX QX 

^l '^l + l 



CV QV 



A5r5f+i, 



-1 < t < 1 
(2.2) 

The magnitude of polymerization, periodicity, anisotropy 
parameter, magnetic field, the electronic g-factor and 
Bohr magneton are represented by t, p, A, il, g and /xb, 
respectively. In the following, we take the unit g^B = 1- 
For t = 0, the low energy sector of this model with 
arbitrary magnetization can be expressed by a Gaus- 
sian model using the Jordan- Wigner transformation and 
bosonization technique as, 



1 



2tt 



vs , 



nG = ^ dxvsKinUr + -^{^) 



K^dx' 



(2.3) 



where (j) is the boson operator restricted to the range 
< (/) < v^TT and n is the momentum density conjugate 
to which satisfies [(/'(a;), n(a;')] = i5{x — x'). The Lut- 
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tinger liquid parameter and velocity of a spin wave are 
represented by K and ws, respectively. This model can 
be described by the conformal field theory with confor- 
mal charge c = 1 . 

If t is finite but small, the low energy sector of our 
model (2.1) with m^ = 1/2 — 1/p reduces to the sine- 
Gordon modcljiJ^ which is given by 



27ra2 



dx 



.^^^(-nf + i(ai) 



/ dx cos v2(/), 



(2.4) 



where the a is the lattice constant, the cos\/2</> term 
generates the energv—gap and its coefficient yi is given 
by 2/1 = 2Trt/vs oc tS3 For small t, yi is small, so that 
this term is regarded as a perturbation in calculating 
the low energy properties. From the flow diagram of the 
lowest order reniarmalization group equation of the sine- 
Gordon modeljIliP we know that there are three critical 
lines which belong to two different universality classes. 
Two are the plateau-non-plateau transition lines K — A = 
2 I 2/1 I . These are BKT transitions. Another is the 
plateau-plateau transition at i = for 1 < K < 4. This 
is a Gaussian transition. The three critical lines merge 
at a multicritical point (t,K) = (0,4). 

2.2 BKT critical point 

It is difficult to estimate precisely the BKT critical 
point from standard finite size analysis of JJie numeri- 
cal calculation data. Nomura and Kitazawao' proposed 
to use the level spectroscopic methodclP with twisted 
boundary condition to overcome this difficulty. This 
method has been successfully applied to the plateau- non- 
plateau transition of the S = 1/2 trimerized XXZ chain 
by Okamoto and Kitazawa.EJ We therefore employ this 
method for the present model. Here we do not explain 
the motivation and background of this method, because 
these are well described in ref. ^^. 

The finite size critical point is determined from the 
crossing point of AiJo,2 and the lower of ^E'^'/^ o defined 

by, 



AEo,. 



i{So(Mp + 2, 0, 1) + EoiMp - 2,0, 1)} 



^o(A^p,0,l). 



and 



AE^ 



/2,0 



= E 



TBC 



(Mp,l)-i?(Afp,0,l) 



(2.5) 



(2.6) 



Ai?i^/2_o = E' ^^ (Mp, -1) - E{M^, 0, 1) (2.7) 

where Eo{M^,k,P) is the lowest energy under periodic 
boundary condition with magnetization M^(=^ Nm'^ = 



N 



Si), wave number k and parity P. The magnetiza- 

1=1 

tion on the plateau is denoted by Mp{= pL{l/2 — 1/p). 

The energy E'^^'-' [M^ , P) is the lowest energy with the 

twisted boundary conditions with magnetization IvP and 

parity P. 

To confirm the reliability of this method, we have also 



checked that the following average 

^1/2, o(-^) '^ 3a;^y2(^) 



(2.8) 



is close to 0.5 at the critical pointstll' where x^J'fj o(^) ^-"^^ 
the scaling dimensions corresponding to £'1)2 n defined 

by 



"1/2,0 



2^^^V2.o- 



(2.9) 



(2.10) 



Here Vg is the spin wave velocity given by 

1;,= lim ^[Em,m{L) - Em,{L)1 

where Em .ki (L) is the energy of the excited state with 
wave number fci = ^ and M^ — Mp. The results are 
shown in Fig. 1 for p = A which confirms that x = 0.5 
holds with good accuracy. 
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Fig. 1. Extrapolated values of the averaged scaling dimension x 
on the critical points for p = 4. The solid line is x = 0.5 



Figure 2 shows the behavior of Ai?o,2 and AEy^ ^ 



for 



L = 6, A = —0.9 and p — 4. From the crossing point, we 
obtain tc{L = 6) = 0.4534. The BKT transition point 
for the infinite system can be obtained by extrapolating 
from L = 2, 4 and 6 to L ^ oo as tc = 0.4537 assuming 
the extrapolation formula tc{L) 
in Fig. 3 for p = 4 and A = 
value is represented by x . 



-0.9. 



■^ as shown 
The extrapolated 



Using the same method, we determine the phase 
boundaries for p = 4, 5, 6, 7 and 8. The phase diagrams 
are shown in Fig. 4. There are two BKT transition 
lines which merge on the multicritical point A — Ac and 
t — 0. The line t = 0, A > Ac is Gaussian corresponding 
to the uniform XXZ chain which has no plateau for any 
value of magnetization. 

From the Fig. 4, we find that the non-plateau re- 
gion becomes narrow and the multicritical point Ac ap- 
praoches —1 with increasing periodicity p. These numer- 
ical results of multicritical points can be checked by the 
Bethe ansatz method as explained in the next section. 
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sented by • and o, respectively, for p 



§3. Multicritical points determined by Bethe 
ansatz method 

By the numerical methods, the p-dependence of the 
multicritical points can not be calculated for large values 
of p. However, the multicritical point is located on the 
line t — which corresponds to a uniform XXZ chain. 
We can therefore make use of the Bethe ansatz solu- 
tion Jiihich is available for arbitrary values of magnetiza- 
tionBB) 

Within the XY-hke region -l<A<l,Ais 
parametrized as A = cos6'. The magnetization m^ and 
Luttinger liquid parameter K are determined in the fol- 
lowing way. First, introduce the function a{ri) for the 
density of particles satisfying the equations, 

a{7^)^^{g{v)- I A{ri^rj')a{ri')d7j'], (3.1) 
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Fig. 3. The extrapolation of tc for p = 4, A = —0.9. 




Fig. 4. Phase boundaries on the A — t-plane for p = 4, 5, 6, 7 and 
8. The solid lines are the BKT plateau- non-plateau critical lines 
and the dotted line is a Gaussian plateau-plateau critical line. 
The lines are visual aids. 



a{r])dT]^-{l~2m'), 



(3.2) 



where the kernel A(ri) and the inhomogeneous term g{ri) 
are given by 

tan 6* 



Aiv) = 



tan^ 9 cosh^ f + sinh^ 2 



giv) 



cosh" ^ 



cot ■ 

I- cot 



2 I sinh 



2 2 
2 



(3.3) 



(3.4) 



The real parameter A > in eq. (3.1) and ( |3.2[ ) describes 
the values of the spectral parameter rj at the Fermi sur- 
face. Solving eq. (3.1), the magnetization m^ can be 
calculated using (|3.2|) for a given value of A. 

On the other hand, the Luttinger liquid parameter K 
is determined by the solution of the integral equation for 
the dressed charge function ^(77). 



av) 



1 

2^ 



A{v-V')av')dv', (3.5) 



as 



K = 2e{A), 



(3.6) 



At the multicritical point, the Luttinger liquid param- 
eter K is 4. For the highest plateaus, the magnetization 
m^ is 1/2 — 1/p. Thus, the anisotropy parameter at the 
multicritical point Ac can be determined from eqs. ( [3.l| ), 
(|3^), ( |3l| ) and (|^. The 1/p-dependence of the mul- 
ticritical points Ac are shown in Fig. 5. The results of 
the Bethe ansatz method and the numerical method are 
represented by o and x , respectively. The results of nu- 
merical calculation coincide well with the Bethe ansatz 
results for 4 < p < 8. 



Investigation of the behavior of Ac for large p is now 
done in detail. Near A ~ — 1, define an infinitesimal 
deviation of A from —1 by i5A = A -I- 1 = cos6' -I- 1 ~ 
(tt - 61)2/2 and (50 = TT - = ^2^. Since A is smah, 
T] and ij' are also small. Therefore, eqs. (3.1) and (3.5) 
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Fig. 5. The 1/p-dependence of the multicritical critical point Ac. 
The results of the Bethe ansatz and the numerical method are 
represented by o and X, respectively. The solid line is the ana- 
lytical result for A — > — 1. The dashed line is a visual aid. 



can be transformed as 
6Q 1 

47r 1-K 



a{x) 



4cr(a;') 



/^'-4+(x-a:')2' 



^(^)^^ + i/ 



dx ^ , 
p 4+ (a; 



4C(x') 



x'Y 



(3.7) 



(3.8) 



where x and /i are defined by a; = r\lbd and /i = K/59. 
Comparing the inhomogeneous terms of eqs. (3.7) and 
(U), we get a{x) ^ ^i{x). From eq. (^j, we find 



l-2m^ 



259 I dxa{x) 



27r 



dxS^ix) 



1 



TT 



dx^{x). 



(3.9) 



By solving eq. (3.S) numerically and adjusting the value 
of /i to satisfy iiT = 4, we can determine the value of /i at 
the multicritical point. Using the thus obtained /i and 
^(z), the relation between Ac and m^ can be calculated 
by (3.9) as 1 - 2m^ = 1.138352(1 + A^) (m^ = 1/2 - 
1/p). The resulting relation between Ac and p is plotted 
in Fig. 5 as a solid line. This result shows that the 
non-plateau region becomes narrow with increasing p but 
does not vanish for finite p. It might appear strange 
that the plateau region is enhanced in the large p limit, 
because our model tends to the uniform XXZ chain in 
this limit. This is not surprising, however, because the 
magnetization at the plateau m^ = 1/2 — 1/p tends to 
the saturation magnetization m^ = 1/2 in the large p 
limit. The highest plateau state just tends to the fully 
magnetized state in this limit. 

§4. Summary and Discussion 

The highest magnetization plateau state at m^ = 
1/2 — 1/p in S* = 1/2 p-merized XXZ chain is inves- 
tigated by exact diagonalization of finite size systems 
and the Bethe ansatz method for general values of p. 
The BKT transition points are determined precisely by 
the level sr)cc±|macppy method with twisted boundary 
condition.tj'ElUlScj-' In the phase diagram, we find two 
plateau- non-plateau BKT transition lines, one plateau- 



plateau Gaussian line and one multicritical point for each 
p. The non-plateau region becomes narrower with in- 
creasing periodicity. 

The multicritical points are calculated by the Bethe 
ansatz method and the results coincide well with the nu- 
merical results for 4 < p < 8. For large p, it is explicitly 
shown that the distance between the multicritical point 
Ac and the ferromagnetic point A = — 1 is proportional 
to 1/p. This implies that the non-plateaus region does 
not vanish for finite p. As the periodicity p becomes 
longer, the highest plateau tends to the fully magnetized 
state. 

The tendency that the plateau region is enhanced with 
the periodicity p is also found in the antiferromagnetic- 
(ferromagnetic)„ polymerized S = 1/2 chains.O From 
these examples, we may speculate that the plateau is 
stabilized as the size of the unit cell increases in gen- 
eral. It would be interesting to investigate a wider vari- 
ety of models with long spatial periodicities to confirm 
this speculation. 

In this work, we confined ourselves to the S = 1/2 p- 
merized XXZ chains with identical strength of modula- 
tion t for the XY and Ising components of the exchange 
coupling. If mapped onto the spinless fermion chains, 
the former corresponds to the transfer modulation and 
the latter, the nearest-neighbour interaction modulation. 
These two are, in principle, independent quantities. In 
fermionic language, the plateau-non-plateau transition 
corresponds to the metal-insulator transition. The non- 
plateau state corrsponds to the state in which the band 
gap is destroyed by the attractive interaction. In this 
context, it must be important to investigate the effect 
of transfer modulation and interaction modulation sepa- 
rately. This is left for future studies. 

The numerical calculation was performed using the 
program package TITPACK version 2 coded by H. Nishi- 
mori on HITAC S820 and SR2201 at the Information 
Processing Center of Saitama University and HITAC 
SR8000 at the Supercomputer Center, Institute for Solid 
State Physics, the University of Tokyo. This work is sup- 
ported by a research grant from the Natural Science and 
Engineering Research Council of Canada (NSERC) and 
a Grant-in- Aid for Scientific Research from the Ministry 
of Education, Science, Sports and Culture of Japan. 
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